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, Abstract. The general line element corresponding to the family of algebraically 

general, gravito-electric, expanding, rotating dust models with one functionally 
independent zero-order Riemann invariant is constructed. The isometry group is 
O ! at most one-dimensional but generically trivial. It is shown that the asymmetric 

solutions with constant ratio of energy density and vorticity amplitude provide first 



, examples of Petrov type I space-times for which the Karlhede classification requires 

I i ' the computation of the third covariant derivative of the Riemann tensor. 

> 

\0 ' PACS numbers: 04.20.Jb 

! 1. Introduction 



Rotating dust solutions in general relativity may serve to describe phenomena on a 



O ■ galactic scale. The metric gab obeys the field equation 

; Rab-^Rgab + ^gab = fJ'UaUb (l) 

c3 with, as usual, Rab the Ricci tensor, R the Ricci scalar, A the cosmological constant, 

■u" the dust 4-velocity field and /i the energy density. For a space-time filled with 
rotating dust the fiuid fiow is non-accelerating and the remaining kinematic variables 
are the expansion scalar 9 = u'^-a, shear tensor aab = U(^a;b) — | hab and vorticity vector 
a;" = |e"'"=U[b;c] 7^ of the fiuid, with hab = Qab + UaUb and eabc = Vabcdu"^ the spatial 
permutation tensor. 

Important classes of rotating dust models have been found by assuming some 
kind of symmetry, or are algebraically special. Respective examples are Winicour's 
classification [1] of stationary axisymmetric models satisfying the circularity condition 
(see [2] for examples and further discussion), and the general rotating dust solution 
admitting time-like conformally fiat hypersurfaces with zero extrinsic and constant 
intrinsic curvature as found by Stephani [3] and later generalized by Barnes for non-zero 
A [Ij, which depends on seven free functions of one coordinate and which turns out to 



I without non-trivial isometries 
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be of Petrov type D. A final and famous example is the homogeneous Petrov type 
D Godel universe [5], which can be interpreted as a rotating dust space-time with a 
negative cosmological constant. 

There seems to be, however, a lack of algebraically general, asymmetric solutions 
In a search for such, the class A of Petrov type /, gravito- electric rotating dust models, 

1. e. for which the Weyl tensor wrt observers comoving with the dust is purely electric, 

Hab = ■^^acdC'^'be = 0, Eab = CacbdU^W^ 7^ 0, (2) 

has been scrutinized [6j. As a main result, it has been proved that the vorticity uj"" of 
u"" must be a geodesic eigenvector of E^-fj. Moreover, the corresponding eigenvalue is 
linearly related to the energy density /i, which cannot be constant, as this would lead 
to a set of Petrov type D solutions containing the Godel universe. Hence, if we denote 
tm for the number of independent components of the Riemann tensor and its first m 
covariant derivatives wrt the Weyl principal tetrad {E"'b eigenframe), we have either 
to = 1 or to = 2. The to = 1 subclass splits into two separate families: the first consists 
of all non-expanding solutions {6 = 0, for which necessarily A < 0), whilst the second 
forms a particular set, say S, of expanding models (^^ 7^ 0), all having A > 0. 

In this communication I deduce the general line element for the family S directly 
from its invariant description wrt the Weyl principal tetrad. It depends on three free 
functions of one coordinate. The construction shows that a metric in S admits at least 
three, but generically four, functionally independent scalar invariants, and thus either 
possesses a group Gi of isometries or is asymmetric. In a next step, the discussion is 
widened to comprise the non-rotating limit case. Finally, denoting q for the number 
of covariant derivatives of the Riemann tensor required in the Karlhede invariant 
classification algorithm [7], it is shown that space-times belonging to a particular 
subfamily of S have g = 3. To the best of my knowledge, this is the highest value 
from Petrov type I examples analysed so far, the previous one being trivial (g = 1). In 
combination with the remarkable recent result by Milson and Pelavas [8], who exhibit 
a set of Petrov type space-times with the theoretically maximal value q = 7 (which 
turn out to be the unique solutions with this property [9]), this reopens the question 
whether the upper bound g = 5 is sharp in the algebraically general well. 

2. Line element and Karlhede classification 

For a generic member of S, we write B = {do = u, 81,82, 83) for the essentially unique 
Weyl principal tetrad, and (fi*^, fi^, fi'^) for the dual basis of one-forms. According 
to the results of [6] we may arrange the tetrad such that the invariant description wrt 
B reduces to the following: 

(i) Curvature variables (/ G M): 

A = A^ fi = 2l{e-l), (3) 

§ To the best of my knowledge, no such models have been found so far. 
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Hab = 0, Ei2 = -^13 = -E'23 = 0, (4) 

Ell — E22 = A^, E22 — E33 = —W, E33 — Ell =l{6 — I). (5) 

(ii) Commutator relations: 

[di, 9o] = [^1, 82] = 0, [9o, ^2] = -/92, (6) 

[^0, d,] = -2uod2 -{e- 1)83, [82, ds] = -2ujdo, (7) 

[di,d,] = -[3d,. (8) 

(iii) Ricci and Bianchi equations: 

doiu = —6uj, diUJ = —I3uj, d2UJ = 0, (9) 

doe = -e{e-i), did = ~f3{e -i), 826 = 0, (10) 

^0/5 = di(3 = -(3^ -W, d2f3 = 0. (11) 



Note from (l3H5l) that there is only one algebraically independent Riemann curvature 
component, and that I ^ 0, 6 ^ and 6 ^ I hj the Petrov type I assumption. The 
describing variables uj = uJi = (723, 6 and (3 are invariantly defined scalars, as they are 
linearly related to commutator coefficients of the geometrically fixed tetrad B. One 
checks that the formal PDE system (iQl fTTI) is consistent with ([6]). The c^s-derivatives of 
UJ, 9 and 13 being undetermined, there will be three free functions in the general solution, 
which we now construct. 

To start with, we note from ([6]) that do, di and 82 forms a Lie subalgebra, i.e., the 
vector field is normal. Hence functions F and z exist such that VL^ = —dz/F. F is 
related to the commutator coefficients 'j^jk via Cartan's first structure equation for fl^, 
dQ^ = —^'j^jkQ^ A Q'', which is equivalent to diF = ^^^F, i < 3, i.e., 

d^F = -{9 -l)F, diF = -(3F, ^2^ = 0. (12) 

Conversely, for any solution F of (1121) it follows that a function z exists such that 
= —dz/F. Fixing one such solution F and defining, for any given scalar invariant 
S, the invariant Sp by 

d^S = FSf, (13) 

the commutation relations [di, d^lS = 7-^,3 djS read 

2 

FdiSp = d^diS + I'is djS, I < 3, (14) 

and the differential d^* can be expanded as 

2 

dS = YdiSn' - Spdz. (15) 

For suitable F, the choice of which we postpone at the moment, the function z 
will serve as one of the coordinates. One observes from (EHH]) that the null vector 
fields 82 ± i9o are also normal and that the normalized vorticity vector di is normal 
and geodesic. However, we will not introduce according coordinates; the deduction of 
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the line element and the further discussion are much more elegant and clear when we 
proceed by constructing the three remaining coordinates directly from the invariants 
and their c^a-derivatives, hereby exploiting the expansion ( ITSj) . 

Firstly, we see from ([9]) and (fTOj) that di{uj/ {6 — I)) = 0, i < 3; it then follows from 
( TTSll that the ratio S = uj/{9 — I) is a function of z, which we henceforth denote by 
1/ f^{z) or, in view of the non-rotating limit to be considered later, by g^i^z). Next, we 
define the invariants t and x by 

p , e 

tanx = — — , cosxe = -. (16) 

— I 

This is equivalent to writing (3 and ^ in a different way, 

/sinxe* ^ /cosxe* 
cos xe^ — 1 cos xe^ — 1 

so that the vorticity amplitude and matter density read 

9-1 1 2P 

Jal^j /3(,-2j(,cosxe* — Ij cosxe' — 1 

Equations (fTOjl and (fTTj) imply the simpler derivatives 

dot = l, dit = 0, d2t = 0, (19) 
^ox = 0, dix = I, d2X = 0. (20) 

We now use ffT^ . Putting 5* = x we get diXp = 0, z < 3, whence xp = fi{z) by ffT51) 
with S = xp. Putting S" = t and rewriting tp = —2y we obtain 

ao|/ = 0, d,y = 0, d2y = y. (21) 

This is the point where we make a convenient choice for F. Looking at (I9j), ( fT2i) and 
(fT9|) we see that we may pick F = cue*, such that d2y = le~^ by (fT2l) . Finally, putting 
S* = ?/ in ([111) we find 

9oi/F = -2/e-2*, = 0, 92?/^ = 2yle-\ (22) 

i.e. = diiy'^ + e"^*), z < 3, such that from (|T5l) with S = yp — {y^ + e~^*) we derive 
l/F = /2(^)+Z/' + e-2*. 

Assembling the above pieces one gets 

dt = m^ + 2ydz, dx = m^ - fi{z)dz, 
dy = le-'n^-{f2{z)+y^ + e~^')dz, = fs{z){e-' - cosx)dz. ^ ' 

Hence, on using z and the scalar invariants t, x, y as coordinates, we infer immediately 
from ( l23i) that the line element is given by 

XMs^ = -{dt- 2ydzf + {dx + fi{z)dzY + (cos x - e~'ff3{zfdz^ 

+ e'^ [dy + {f2{z) + + e~'') dz]\ (24) 

where the invariant scalar fields f\{z), f2{z) and fsiz) are arbitrary functions of their 
argument. Notice that A = a/A plays the role of a constant scaling factor. The variable 
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t is time-like and x and y are space-like, whereas z is spacelike, null or timelike in the 
region where 

9.. = fl + e'\h + e-^* + y'f + (cosa: - e'^ fl - V (25) 

is greater than, equal to or smaller than zero, respectively. Switching to the null 
coordinates u± and the space-like coordinate ^ defined by 

u± = y±e-\ e = x-0(z), ^(2;) = -/i(2;), (26) 
the metric becomes 

AMs^ = (^2 + ^0^{^2 _ ^0) _^ ^^1)2 ^ ^^3)2^ (27) 

f]2^f]0_ 1 [d«^ + (/2(;2)+4)dz], (28) 

VLj^ — VL — 

n' = di, n^ = f,{z){co^{i + <P{z)) + '^^^^dz. (29) 

(j28l) makes the normality of =F apparent, as du± + {f2{z) + u'^)dz are one-forms 
on 2-spaces; at the same time it reveals the difficulty when we would have started with 
coordinates v± and scalar fields G±(f+, f-, ^, z) for which fi^ =f = G±dv±. 

One observes from f[T7|) and ffTS]) that yU > 0, and at the same time > 0, at 
space-time points with cosx > e~*, while the boundary cosx = e~* forms a two-surface 
of curvature singularities. It is an open question whether ([21]) may generate physically 
plausible rotating dust models. 

The discussion of possible Killing vector fields becomes trivial in the above invari- 
ant approach: when at least one of the functions fi{z), f2{z) or fz{z) is non-constant it 
can be taken as a fourth invariantly defined coordinate replacing 2, and there is no sym- 
metry; when all of them are constant there is a group Gi of isometries generated by d/dz. 



By putting u equal to zero in fl MTTj) one arrives at the corresponding irrotational 
dust family , which was integrated in [10] (metric (45)). More generally, we may 
consider the class S = SUS^ . Then the above integration procedure can be applied up 
to ( |2T]) . and we make the alternative choice F = [9 — l)e^ . Taking further advantage of 
the normality of e.g. 82 — do and solving Cartan's first structure equations, one readily 
obtains the line element 

Pds^ = - (dt - 2ydzf + {dx + fi{z)dzf + (cosx - e-'fdz^ 

+ {e^-^dri + dt-2y dzf, (30) 



where now 



y = y{t, V, ^) = ^) " ^ ^33iz) (31) 



and where K = K{r], z) is a solution of 



^''^'''^^|^(^,^)=2^73(^). (32) 
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The non-rotating subcase corresponds to gs^z) = 0. Looking at (ETj) . (15T]) and 
fl32|) . and after possibly redefining rj, we have that K = K{z) is a primitive function of 
2y = 2y{z). In terms of coordinates z, t = t — K{z), ( = ri — e~* and ^, cf. ( l26l) . the 
non-rotating solutions take the form 

AMs^ = - dr^ + de' + e^^C' + [e-^(^)e-" - cos{^ + (j){z))]^dz^ . (33) 

There is at least a Gi isometry group generated by 9^ ~ e'^d2- As the scalar invariants t, 
X, y and /i must be invariant under any isometry, one immediately deduces that there 
is a group G2 of motions if and only if both y and /i are constant, i.e., if and only if 
the line element can be transformed to fl33l) with K{z) = yz and (f){z) = fiz, where y 
and /i are constants; the additional Killing vector field is then given by 

2ydr - -dz-^ 2ydo - fidi + (e"* - cos x)d3. (34) 
This corrects an inaccuracy in the symmetry discussion of 5° in [10] . 

When working wrt the Weyl principal tetrad i3 of a Petrov type / space-time, the 
stop value q within the Karlhede invariant classification algorithm is defined by 

<m < q, tq = (35) 

where one formally puts t_i = —1. The final value tg_i is the total number of 
functionally independent scalar invariants and equals 4 — r, where r stands for the 
dimension of the isometry group. For members of S one has t^ = 1 and tg_i = 2, 3 or 
4, such that 2 < g < 4 a priori. A direct calculation shows that the components of 
the first covariant derivative of the Riemann tensor wrt B are rational functions of the 
invariants 

cosa;e*, sinxe*, 2 ?/ cos a; + sina; /i(2;), ga^z). (36) 
The following possibilities thus arise (/' = means that / is constant): 

• for metrics in {gsiz) = 0, y = y{z)): 

(i) fi = y' = 0. There is a group G2 of isometrics and ti = 2, such that t2 = 2 
and q = 2. 

(ii) /( 7^ or 7^ 0. There is a group Gi of isometrics and ti = 3, such that t2 = 3 
and q = 2. 

• for metrics in S (gsiz) = 1/ /^{z) 0, y functionally independent of t, x and z): 

(i) f[ = f2 = fs = 0. There is a group Gi of isometrics and ti = 3, such that 
^2 = 3 and q = 2. 

(ii) /g = 0, but /( 7^ or /2 7^ 0. There is no symmetry and ti = 3, such that 
h = h = 4: and q = 3. 

(iii) 7^ 0. There is no symmetry and ti = 4, such that ^2 = 4 and q = 2. 

We conclude that q = 2 for all members of S, except for the asymmetric (whence 
rotating) models for which the ratio 21 = /i/tu of energy density and vorticity 
amplitude is constant. This special class of models, depending on two free functions 
of one coordinate and one constant parameter, provides a first example of algebraically 
general space-times for which g = 3. 
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